Acquire information about neutrino parameters by detecting supernova neutrinos 



Ming- Yang Huang 1 *, Xin-Hcng Guo 1 ''', and Bing-Lin Young 2, 3 * 

1 College of Nuclear Science and Technology, Beijing Normal University, Beijing 100875, China 

2 Department of Physics and Astronomy, Iowa State University, Ames, Iowa 5001, USA 

3 Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100190, China 

We consider the supernova shock effects, the Mikheyev-Smirnov-Wolfenstein (MSW) effects, the 
collective effects, and the Earth matter effects in the detection of type II supernova neutrinos on 
the Earth. It is found that the event number of supernova neutrinos depends on the neutrino 
mass hierarchy, the neutrino mixing angle 613, and neutrino masses. Therefore, we propose possible 
methods to identify the mass hierarchy and acquire information about #13 and neutrino masses by 
detecting supernova neutrinos. We apply these methods to some current neutrino experiments. 
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I. INTRODUCTION 



Supernovas (SNs) are extremely powerful explosions in the universe which terminate the life of some stars 
They make the catastrophic end of stars more massive than 8 solar masses, leaving behind compact remnants such as 
neutron stars or black holes. The SN explosion, such as SN1987A Q, is one of the most spectacular comic events and 
hH ■ the source of new physics ideas Q . Detecting SN neutrinos on the Earth has been a subject of active investigation in 
astroparticle physics [l| Hfl- S ■ We may obtain some information about the explosion mechanism of SN Q and the 
Q_C parameters of neutrinos 0] [H by detecting SN neutrinos on the Earth. Several neutrino experiments, including 
the Daya Bay reactor neutrino underground laboratory Q , can be used to detect neutrino events from a SN explosion 
and serve as a SN Early Warning System [h|. 

Recently, several methods to obtain information on neutrino parameters by the study of SN neutrinos have been 
proposed. In Ref. 0], Dasgupta et al. proposed a method to identify the neutrino mass hierarchy for very small #13 
l— — '■ (sin 2 #13 < 10 -5 ) through the Earth matter effects. In Ref. 0], Skadhauge et al. proposed methods to identify the 
I neutrino mass hierarchy for larger #13 (sin 2 6*13 > 10 -4 ) and obtain information about 6*13 by SN neutrinos. In Ref. 0], 
■ using the event number of the delayed SN neutrinos, a method to measure neutrino masses was proposed by Beacom 
£■ — ' et al.. In these works, the SN shock effects and the collective effects were not taken into account. Furthermore, 
0^ ■ uncertainties of some of the parameters appearing in SN neutrino fluxes, such as the luminosities lP a , temperatures 
T a , and pinching parameters rj a , were not considered. 

In Ref. , we studied the Earth matter effects in the detection of SN neutrinos at the Daya Bay experiment and 
CO ' proposed a possible method to acquire information about the neutrino mixing angle #13 smaller than 1.5° which is 
] difficult to access in reactor neutrino experiments. We defined the ratio R as the event number of v e over that of u e . 
Using the relation between the event number of SN neutrinos N and #13, we can obtain information about small #13 
by measuring R. In that paper, a simplified picture for the collective effects was used [l2j and the neutrino flavor 
transformation due to the SN shock effects was not considered. Furthermore, it was assumed that different flavor SN 
neutrinos have the same time- independent luminosity. Hence in that calculation there are several uncertainties which 
we will address in the present work. 

In this paper, we study the collective effects more accurately [Hjj[l4| and consider the neutrino flavor transformation 
arising from the SN shock effects [l|| [l?} . In addition, according to current Monte Carlo simulations on SN neutrino 
spectra (l8j . we allow different flavor neutrinos to have different time-dependent luminosities. Furthermore, in order to 
reduce the uncertainties of unknown parameters in SN neutrino fluxes such as luminosity, temperature, and pinching 
parameter, we will define the ratio of the event number of SN neutrinos for the charged-current reaction detected 
after one second over that detected before one second (this is different from the ratio considered in Ref. [HI). With 
this ratio, we may identify the mass hierarchy, Am 2 ^ > vs Am 2 ^ < 0, and acquire information about #13. This 
method can be applied to neutrino detectors of various types. Examples are the liquid scintillator detector of the Daya 
Bay experiment [9| under construction, the decommissioned heavy water detector of the previous Sudbury Neutrino 



* Email: hmyl9151905@mail.bnu.edu.cn 

T Corresponding author, email: xhguo@bnu.edu.cn 

* Email: young@iastate.edu 



2 



Observatory (SNO) [H, etc. 

The other main purpose of the present work is to propose a method to obtain some information about neutrino 
masses from SN neutrinos. Since neutrinos have small masses, their transmit velocities are smaller than that of light. 
Therefore, a massive neutrino will have an energy and mass dependent delay relative to a massless neutrino after 
traveling over a long distance @. We will define another ratio, for the neutral-current reactions, as the event number 
of the delayed SN neutrinos over the total event number. With the relation between this ratio and neutrino masses, 
we will propose a possible method to acquire information about neutrino masses. This method will also be applied 
to Daya Bay and SNO for examples. 

The paper is organized as follows. In Sec. II, we review the detection of SN neutrinos on the Earth. In Sec. Ill, we 
propose a possible method to identify the mass hierarchy and acquire information about # 13 and apply the method 
to the Daya Bay and SNO experiments. In Sec. IV, a method to acquire information about neutrino masses is given 
and also applied to Daya Bay and SNO. A summary is given in section V. 



II. DETECTION OF SN NEUTRINOS ON THE EARTH 



In the explosion a typical type II SN generates intensive neutrinos which carry about 99% of the total energy and 
the explosion itself consumes about 1% of the total energy. The vast amount of neutrinos are produced in two bursts. 
In the first burst which lasts for only a few milliseconds, electron neutrinos are generated by the electron capture by 
nuclei e~ + N(Z, A) — > N(Z — 1, A) + v e via the inverse beta decay e~ +p — > n + v e which leads to a neutron rich star. 
In the second burst which lasts longer (about 10 seconds), neutrinos of all flavors are produced through the electron- 
positron pair annihilation e~ + e + — > v a + P a , the electron- nucleon bremsstrahlung e + -/V — > e + N + v a + D a , the 
nucleon-nucleon bremsstrahlung N + N — > N + N + is a + D a , the plasmon decay 7 — > v a + v a , and the photoannihilation 
-1 + e ± e± + v a + D a [U. 

We assume a "standard" SN explosion at a distance D — lOkpc from the Earth, releasing a total en ergy Eb = 
3 x 10 53 erg (similar to SN1987A) and the luminosity flux of the SN neutrino a, L a , distributed in time as [2lj [22l|[23j 

L a (t) = L° a e-*'T, (1) 
where L° is the luminosity of the neutrino a at t = and r the luminosity decay time. The following rang e of r was 



obtained by fitting the the experimental data of SN1987A with the exponentially decaying luminosity [2l[[22j: 

T = I.74 - 4.19s. (2) 
The typical relations among L° obtained from numerical simulations are [lH 

^ = (0.5-2), L° Ve =Ll, L° Vw =L° mi (x = M,r). (3) 

There are two parametrizations of SN neutrino fluxes used in general simulations, one used by the Livermore group 
[24| [lU and the other by the Garching group [l8| . In the present work, we use the Livermore parameterization which 
is similar to the Fermi-Dirac distribution. For the SN neutrino of flavor a, the time-integrated neutrino energy spectra 
is given by [2f| 

F " 0) {E) = F^Tl cxp(E/T a - Va ) + l ' (4) 

where E is the neutrino energy, T a the temperature of the neutrino a, rj a the pinching parameter of the spectrum, 
L a the luminosity given by Eq. (JlJ, and F a j is defined by 



where j is an integer. The spectra obtained from numerical simulations can be well fitted by (26| 

T Ve = 3 - 4Afe\/, T Pc = 5 - 6M eV, T^ x = T Vi = 7 - 9M eV, (5) 
rive w 3 - 5, ??Pc Ps 2.0 - 2.5, r)„ x = ??Px w - 2. 

While neutrinos propagate outward to the surface of the SN they could be subjected to the SN shock effects 
HI (H [I!] (H, the MSW effects [11 [U [HI , and the collective effects E3l H3| H^l H3| - Then, 



after travelling the cosmic distance to reach the Earth, they go through a part of the Earth and are subjected to the 
Earth matter effects gj g| . In the following, we will consider the above four types of effects in detail. 
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A. Time dependence of the signal: SN shock effects 

As was pointed out in Ref. [lj| , the shock propagating inside the SN may reach the region of densities relevant for 
neutrino conversion during the period of neutrino emission. It modifies the density profile of the star, thus affecting 
neutrino flavor transitions after the core bounce 1|| 17]. The time-dependent variations of the neutrino potential in 



matter along with the SN shock profile can leave distinctive imprints on the energy and time structure of the neutrino 
signal. The study of such effects may provide us with additional constraints for the determination of the neutrino 
mass hierarchy, the neutrino mixing angle #13, and neutrino masses. 

In Ref. [l7J, the authors introduce a simplified, empirical parametrization of the shock density profile p, which 
reproduces the main features of the graphical profile of Ref. [la ] and is continuous in both t and the supernova radial 
coordinate r. In this way, numerical and analytical calculations of the crossing probability at the high resonance 
region Ph can continuously cover the relevant parameter space. We follow the notation of Ref. [l7j in the discussion 
below. 



For post-bounce time t < Is [44J before the shock effects begin, the SN matter density is 

Po(r) * 10 14 • ^jj-^ g/cm 3 . (6) 
For t Is the shock effects set in and the matter density is given by 

p(r,t) = po(r) ■ < / \ (7) 

[1 {r>r s ), 

where f(r,t) is defined as 

f(r,t) = cxp{ [0.28- 0.69 ln(r s /fcm)][arcsin(l - r/r,)] 1 - 1 }, (8) 
and £ is a typical ratio of the potential across the shock front, 

£ = V+/V. - 10, (9) 
which measures the SN matter potential V(r) drop from 

V+ = lim V(r) (10) 

t— >r B 

to 

V- = lim V(r). (11) 

r—tr£ 

The SN matter potential is related to the SN electron density N e (r) by 

V(r) = V2G F N e (r) = V2G F N A p(r)Y e , (12) 

where Na is the Avogadro's number and Y e is the electron fraction. In the numerical calculations, we assume Y e = 0.5. 
In Eqs. a slightly accelerating shock-front position r s is assumed with the explicit time dependence, 

r s (t) = -4.6 x 10 3 + 1.13 x 10 4 • t + 1 x 10 2 • t 2 , (13) 

where r s is in units of km and t in units of s. 

Using Eqs. (jU)-©, we plot in Fig. 1 the density profile of the SN as a function of the radius at different times. 
These curves reproduce reasonably well the main features of the shock profiles as shown in Ref. [l6| . The straight line 
labelled by t = is the unmodified density profile given by Eq. ^ . Once the shock effects begin, at a given time the 
shock radius r s is determined by Eq. (fT3")) . The density profile is suppressed more and more in the region r < r s as the 
time and therefore r s increase. The density reaches a minimum and then rises relatively quickly to the value £/°o( r s) 
as r — > r~ . In the region r > r a the density profile is the unperturbed original form pq. The density suffers a sizable 
discontinuity at r^. 1 For neutrinos emitted in the region r > r s , due to their much higher velocities, which is the order 



1 It can be shown that f(r s ,t) = 1 and therefore the discontinuity is (£ — l)po(r s ) Kb 9p(r a ). 
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of speed of the light, than that of the shock wave propagation, which is the order of 10 4 km/s, they are not influenced 
by the shock effects. Hence the majority of SN neutrinos which are emitted within a 10 second period will be able to 
semble the whole range of the profile of the SN density under the influence of the shock. An important consequence 
of the density modification by the shock wave is that some of the neutrinos can undergo multiple resonance flavor 
conversion when they traverse the SN. We consider them below. 

Let us first consider the neutrino flavor conversion due to the density discontinuity at the shock front r s . This 
crossing probability of the neutrino flavor transformation which we denote as P s can be determined by the conservation 
of flavor across the matter potential discontinuity. By a simple calculation, we have (3l| 

P s =sm 2 (e+-6-), (14) 

where and #~ are respectively the effective mixing angles in matter immediately before and after the shock front 
position r s . Such angles are defined by 

cos2^ = - Ar^ cos 2ft, -21^* 

^(Am§ x cos26» 13 - 2EV ± ) 2 + (Am^ sin26» 13 ) 2 



B. MSW effects 

The MSW effects [3(| is the adiabatic or partially adiabatic neutrino flavor conversion in media with locality or 
time varying matter density. The main feature is that in a high density medium a flavor is also a mass state which 
propagates without flavor oscillation, until crossing a resonance where there is a finite probability for flavors change, 
i.e., a jump from one mass eigenstate to another. For #13 smaller than 3°, using the Landau-Zener formula [45|, 
one can obtain the crossing probability, Pm, for neutrinos to jump from one mass eigenstate to another at the high 
resonance layer (r = r,) 

P m = exp(-^y), (16) 



with 



LA^sin^2gi3 1 
/ 2E cos 26>i3 I d In N e /dr | res ' 1 ' 



where N e is the electron density 31J and jAm^l ^ |Ato3 2 | ^ 2.4 x 10 3 eV 2 [46]] ■ Similarly, one can calculate the 



expression of the crossing probability Pj, at the low resonance region inside the SN. However, due to the large angle 
solution of the solar neutrino, Pl ~ 0. 

In the presence of a propagating shock wave, the calculation of the global crossing probability, Pjj, is not straight- 
forward. The shock modifies the density profile as discussed in the preceding subsection. Within the shock wave and 
near the shock front the matter density is no longer monotonically decreasing, but rather goes through a minimum 
and a maximum. So the resonance condition can be satisfied possibly as many as three points: 

Am^ ■ cos28 13 = 2y/2G F N e {n)E, (i = 1, 2, 3). (18) 

It is, of course, possible that there are only two or one resonance points, or even none at all. In the case of three 
resonance points, which we denote as rj, i = 1, 2, 3, two of the points lie below the shock front and the third one 
beyond the shock front, i,e., r\ < < r s < r 3 . The crossing probability Pm at r,, i = 1, 2, 3, can be calculated from 
Eq. (fTB")) . Note that although the matter densities are the same at all resonance points as required by Eq. (fl~8|) , their 
crossing probabilities, as given by Eq. (jTTJ) , may be different because the slop of the matter density may be different 
at different resonance points. The global crossing probability can be expressed in terms of four crossing probabilities 
Phi, Ph2, P s , and P H3 . 

To obtain the expressions of the global crossing probabilities, we first define four densities: p res is the resonance 
density which satisfies Eq. (|18|) . p+ the density at r = r s corresponding to the matter potential V+, p- the density 
at r = r s corresponding to the matter potential V— , and pb the density at the bottom of the camber which is the 
minimum of the density profile below the shock front. From Fig. 1, it can be seen that 

P+ > P-, P+ > Pb- 

In the pre-bounce time period t < Is, the resonance condition occurs at only one point and there are no neutrino 
shock effects. Then. 



Ph = Pm- 



(19) 
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In the post-shock period t > Is we have four different situations: 

(1) For p_ < p res < pt,, the resonance condition doesn't occur and then 

P H = P S . (20) 

(2) For p res > p + , or p res = pb and p res > p_, or p res < pb and p res < p_, the resonance condition occurs at only 
one point, r — 7*1 or ?* = r 3 , and then 



Pjj = Pfffe +P.s- 2P Hk P s + 2^P Hk P s (l - P Hk ){\ - P s ) cos</> fcs , (fc = 1,3). (21) 

(3) For pb < Pres < P+ and p res > or p res = and p res < p_, the resonance condition occur at two points, 
r = r\ , r 2 or r = ri , r 3 , and then 

Pff = {Phi+Phi+Ps)-2(PhiPhi+PhiPs + PhiPs)+4PhiPhiP s 

+2(1 - 2P s ) v /P m P Hi (l - -Pm)(l - Pm) cos0 i; (22) 

+2(1 - 2P H iWPhiPs(1 - - P«) cos<fe 

+2(1 - 2P m )y/P H1 P a (l - P H1 )(1 - P s ) cos0 ls (i = 2,3). 

(4) For pb < p res < p-, the resonance condition occur at three points, r = ri,r 2 ,r3, and then 

Ph = (Phi + Phi + P s + Pjjs) - 2(P m P H 2 + fti^ + ftift + Pff 2 P ff3 + Pff 2 ft + Pffsft) 
+4(PffiPff 2 Pff3 + PhiPh2P s + PhiPhsPs + Ph2P H3 Ps) - 8P H iPh2Ph 3 Ps 
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>s3- 



SVPhiPh2(1 - Phi)(1 - Ph2) cos0 12V /P H3 P s (l - P H3 )(l - P s )cos, 
In Eqs. (|21ti23[) . fa is defined as 

<hi ~ y ' da;^^/[Am§ 1 cos2ei3 - 2£F(r)] 2 + (Am^ sin 2# 13 ) 2 , (24) 

with V(r) being the SN matter potential. A summary for neutrino flavor conversions due to the neutrino shock effects 
and the MSW effects in various density regions is given in Table I. Due to the SN matter effects and the SN shock 
effects, the expressions of P H , Eqs. (fl9 ]l -([2"3" |) . are different from Eq. (25) in Ref. [13] and Eq. (34) in Ref. [13]. 

Because of the various effects they have been subjected to, the SN neutrinos reaching the Earth have a very rich 
structures to be explored in their spectra in time, energy, and the mixing angle O13. In Fig. 2, we plot the crossing 
probability Ph as a function of the neutrino energy E at t = 4s for two very small mixing angles sin 2 613 = 10~ 5 
and 1CT 4 . It is seen that there are large variations in Pjj as the neutrino energy varies. These forms of behavior are 



different from those given in Refs. [llj and [17]]. In Fig. 3, we plot Pjj as a function of t for sin #i 3 = 10 at a 
fixed E for three different values of E. It can be seen that the shapes of the curves depend strongly on E and there 
are many oscillations in Pjj in an interval around t = 5s. The interval enlarges and moves to higher t value as E 
increases. This is different from that given in Ref. [13] • In Fig. 4, we plot Pjj as a function of O13 for different values 
of E and t. It is found that at some values of #i 3 , which depends on the energy of SN neutrinos, Ph begins to change 
very rapidly with 613 at t = 4s, and the behavior of the Ph dependence on 613 is different from that given in Ref. 
fill ]. These differences are attributable to the different SN matter effects when the shock effects is included in the 
consideration. The relation between Ph and 613 will provide a way to obtain information on #i 3 and the study of the 
shock wave effects. 



C. Neutrino-neutrino interactions: collective effects 



Because of the neutrino-neutrino forward scattering or neutrino self-interaction, neutrinos can experience collective 
flavor transformation in the SN near the neutrinosphere where the neutrino number densities are very large [35| . This 
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phenomenon is different from the MSW effects in that the flavor evolution histories of neutrinos in the collective 
oscillations are coupled together and must be solved simultaneously. 

Recentl y, a sig nificant amount of studies on collective effects have been made by a number of authors, e.g., Duan 
et al. [33j [34] 35] |3(| , Dasgupta et al. @ [l2j , Hannestad et al. [13] , Raffelt et al. [iH , and Fogli et al. ■ Defining 
Pyv(Pvv) as the probability that the (anti-)neutrino v(y) remains as v{y) after the collective effects, Dasgupta et al. 
[a] introduced a much simplified picture to explain the characteristics of the collective effects which we followed in 
our earlier paper (Tlj . i.e. P vv = P vv = 1 in the case of normal mass hierarchy; and P vv = f , 



p — 



(E < E c ), 
(E > E c ), 



in the case of inverted mass hierarchy, where Ec is a critical energy across which a stepwise flavor conversion of v t 



develops. However, from the numerical simulations of the collective effects obtained by Duan et al. |f4j, it is found 
that P VV (P VV ) depends on the neutrino (anti-neutrino) energy E, the SN radial coordinate r, the SN neutrino emission 
angle, and it is very difficult to obtain an analytical expression for P VV (P VV ). It can be seen from Movies 5 and 6 
in Ref. [TH 2 that when the SN radial coordinate r increases, there are oscillations between v e ( v e ) and v T (v T ) due 
to the collective effects. When r > 220fcm, which is far away from the neutrinosphere, the collective effects for the 
neutrino amount to a swap between the v e and v T fluxes, irrespective of the mass hierarchy being normal or inverted. 
For the anti-neutrino, in the case of normal mass hierarchy, the fluxes of v e and v r are nearly unchanged and there is 
no swap between them. In the case of inverted mass hierarchy, the behaviors of the v e and v T fluxes are very complex 
and, for simplicity, they are generally approximated flU] (35| as that there is no swap between them. In this way, we 
apply the following new simplified picture to describe the characteristics of the collective effects: 

Mi S>K 

for neutrinos and 

Pw = l, (26) 

for antineutrinos. It will be seen from our numerical calculations that the results obtained in Sees. Ill and IV hardly 
change for Ec varying between 6MeV and lOMeV. Therefore, similar to Ref. [f3|, we take Ec = lOMeV. 



D. Earth matter effects 



Neutrinos from a SN may travel through a significant portion of the Earth before reaching the detector and are 
subjected to the matter effects [4(j. Therefore, we need to consider the Earth matter effects in studying SN neutrinos. 
Suppose a neutrino reaches the detector at the incident angle 9 (see Fig. 5). Then the distance that the neutrino 
travels through the Earth to the detector L and the distance of the neutrino to the center of the Earth x are given by 



L = {-R E + h) cosO + yjR 2 E - (Re - h) 2 sin 2 0, 

x = y/{-R E + h) 2 + {L- x) 2 + 2(R E - h){L - x) cosfl, 

where h is the depth of the detector in the Earth, x the distance that the neutrino travels into the Earth, and Re the 
radius of the Earth. 

Let Pi e be the probability that a neutrino mass eigenstate Vi enters the surface of the Earth and arrives at the 
detector as an electron neutrino v e , we have [42j 

P 2e = sin 2 0u + I sin 2 29 12 f ' dxV(x) sin 0™ , (27) 

1 Jxo 



2 Movies 5 and 6 which describe the full simulation of the collectives effects 1 1 -ll are available at 
|http://iopscience.iop.org/1749-4699/l/l/015007/media"| 
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where V(x) is the potential v e experiences due to the matter density p{x) inside the Earth |42] [47] 

V{x) = V2G F N AP {x)Y e , (28) 

and <j>™_+ b i s defined as 



b 

dxA m (x) 



where 

2E 

with 6»i 2 - 34.5°, Amlj = 7.7 x l(r 5 eV 2 , e(x) = 2.EV(:r)/Aml 1 g 



Am 2 / 

A m (x) = -^i^/(cos20 12 -e(a;)) 2 + sin 2 20 12 , (29) 



E. Summary of all matter effects 

When all the effects, including the shock, MSW, collective, and the Earth matter, are taken into account, the SN 
neutrino fluxes at the detector can be written as 

Fg = P fW + (1-P)f£\ 
2F» = ( 1 - p )F^ + (l+p)F^, 

2F° = (l-p)F^ + (l+p)F^ 1 (30) 
where the survival probabilities p and p are given by 

P = P2e[P H Pvv + {l- P H ){1- Pvv% 

p = (l-P 2e )P vu , (31) 

for the normal mass hierarchy and 

P = P2eP V v, 

P = (1 - P 2 e)[P H P,u + (1 - P H )(1 - Puu)], (32) 

for the inverted mass hierarchy. 

Therefore, the event numbers N{i) of SN neutrinos in the various reaction channels "i" can be calculated by 
integrating over the neutrino energy, the product of the target number Nt, the cross section of the given channel 
<r(i), and the neutrino flux function at the detector F® , 

N(i) = N T JdE-a(i)-^-F°, (33) 

where a stands for the neutrino or antineutrino of a given flavor, and D is the distance between the SN and the Earth. 
In the following section, we propose a method to identify the mass hierarchy and acquire information about #13 by 
detecting SN neutrinos. 

III. IDENTIFY THE MASS HIERARCHY AND ACQUIRE INFORMATION ABOUT 6 1Z BY 

DETECTING SN NEUTRINOS 

In the above section, we have discussed the SN shock effects, the MSW effects, the collective effects, and the Earth 
matter effects in detail for the detection of SN neutrinos. From Eqs. (fT4"]) - (fT7]) and (fH)]) - (f3l?]) . it can be seen that 
the event number of SN neutrinos depends on the mixing angle #13 for different mass hierarchies. Therefore, we may 



identify the mass hierarchy and acquire information about # 13 by detecting SN neutrinos. In Rcf. 11], we defined a 
ratio of the event number of v e over that of v e and studied the method to acquire information about # 13 . However, 
in that paper, we ignored the SN shock effects and only a very simple model for the collective effects was considered. 
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In addition, with the naive luminosity of SN neutrinos employed there, the uncertainties in #13 will be very large. 
In the present work, we take into account the SN shock effects and consider the collective effects more realistically. 
In order to further reduce the influence of the uncertainties in the parameters of the SN neutrino fluxes, we define a 
more appropriate ratio as the event number of SN neutrinos for the charged-current reaction detected after Is over 
that detected before Is. Us ing this new ratio, we can better acquire information about #13 and identify the neutrino 
mass hierarchy. From Ref. and detailed numerical studies related to the present work, we found that the ratios of 
event numbers of SN neutrinos are nearly independent of the SN neutrino incident angle. Therefore, in the following 
we will only consider the case for the incident angle of 30° . 

There will be eight detectors located at the near and far sites of the Daya Bay experiment Q. The total detector 
mass is about 300 tons and the depth of the detector h ±2 400 m. The Daya Bay Collaboration has decided to use 
Linear Alkyl Benzene (LAB) as the main part of the liquid scintillator. LAB has a chemical composition including 
C and H and the ratio of the number of C and H is about 0.6. Then, the total numbers of target protons, electrons, 
and 12 C are 

= 2.20 x I0 31 , = i- 01 x 1q32 > N t ] = !- 32 x 1q31 - 

In the Daya Bay experiment, there are three reactions which can be used to detect SN neutrinos: the inverse beta 
decay, neutrino-electron reactions, and neutrino-carbon reactions. However, since the event number observed in the 
channel of neutrino-electron reactions at Daya Bay is very small [ll| , these reactions will not be considered in the 
present work. Therefore, we only consider the inverse beta decay and the neutrino-carbon reactions, their effective 
cross sections are given as follows 48| 49j | : 



(1) The inverse beta decay is 



a{D eP ) = 9.5 x I0- 44 (£(MeF) - I.29) 2 cm 2 , (34) 

where the reaction threshold is Eth = 1.80MeV. 
(2) The neutrino-carbon reactions are 



(a( 12 C(v ei e-) 12 N)) = 1.85 x I0- 43 cm 2 , 

(a( 12 C(D e ,e+) 12 B)) = 1.87 x lO^cro 2 , (35) 



for the charged-current capture, and 



" 43 cm 2 , 



{a{ufC)) = 1.33 x 10 
{(r{v 12 C)) = 6.88 x I0~ 43 cto 2 , (36) 
(^KtOW) 12 <?)) = 3-73 x 10~ 42 cm 2 , 



for the neutral-current capture. The effective cross sections in Eq. (|35j) are valid for SN neutrinos without oscillations. 
When neutrino oscillations are taken into account, the oscillations of higher energy v x into v e result in an increased 
event rate since the expected v e energies in the absence of oscillations are just at or below the charged-current reaction 
threshold. This leads to an increase by a factor of 35 for the cross section (cr( 12 C(^ e , e~) 12 N)) if we average it over 
a v e distribution with T — 8MeV rather than 3.5MeV. Similarly, the cross section (a( 12 C(u e , e + ) 12 B)) is increased 
by a factor of 5. 

In order to reduce the uncertainties of the unknown parameters appearing in the SN neutrino fluxes, including 
those in the luminosities L° a , temperatures T a , and pinching parameters rj a , we define the ratio R,f^ ' as the event 
number of v e (p e ) for the charged-current reaction detected after Is over that detected before Is, 

_ A^ )( ls<t<10s) 
1 " ^)(0 <*<!-) ■ ( } 

With Eqs. ©, O, (fI3 P -(fI7 1) . (fT9 > (|3"3" ]) , and (|37)l. we can see that Rf depends on the mixing angle 6 13 and the 
luminosity decay time r. Form a detailed numerical calculation, we found that Rf which changes with #13 has also a 
significant dependence on r. This can lead to large uncertainties in Rf when r varies in the range given in Eq. ([2"]l. 
Therefore, it will be difficult to identify the mass hierarchy and acquire information on # 13 from Rf when r has a 
significant uncertainty. However, if the value of Rf is determined through SN neutrinos, we can obtain a rather tight 
relation between r and #13. We demonstrate this relation in Fig. 6, for the neutrino-carbon reactions. In Fig. 6 we 
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plot t as a function of #13 for a fixed value of R^, We take Rf = 2.4 and L° = 1 as an example 3 . From this 

figure, we see that the value of r lies in a small range, i.e., 2.68 — 3.00 when (? 13 varies from to 3° for the normal 
mass hierarchy. This is much smaller than the presently allowed range given in Eq. ((2|). It should be remarked that 
"l7( | have typically used the best fit value, i.e., r = 3s [22j in their study of SN neutrinos. 
First consider the channel of the inverse beta decay in which only is relevant. We plot in Fig. 7 as a 
function of #13 for the two mass hierarchies, and for L v jh v , T ai and r\ a taking their allowed extreme values in the 
ranges given in Eqs. ([3]) and (|5|). For the normal mass hierarchy R± , having the value of 2.4, is independent of #13. 
So it is impossible to obtain any information on #13 for the normal mass hierarchy in the inverse beta decay process. 
For the inverted mass hierarchy, for a given value of i Pc /i Px , R^ lies in a band bounded by the curves marked as 
" inverted & max" and " inverted & min" , which are obtained by setting the values of both T a and rj a at their upper 
and lower bounds, respectively. From this plot, it can be seen that if the value of is away from 2.4, that is statically 
significant, the mass hierarchy must be inverted and a band of values of #13 and L Pc /L P;r can be determined. But 
it is difficult to determine #13 from R± alone. It is interesting to observe that, since the value of R± falls bound in 
the limited range of 2.1-2.8, it can serve as a test of the validity of the overall parametrization of the SN explosion 
mechanism. 

For the neutrino-carbon reactions, both R^ and can be used to identify the mass hierarchy and obtain infor- 
mation on #13 . In Fig. 8a, we plot Rf vs #13 similar to Fig. 7. Here for the inverted mass hierarchy, Rf has the 
constant value 2.4 independent of #13. If Rf is larger than 2.4, the mass hierarchy must be normal. It can also be 
seen from Fig. 8a that R± depends on #13 sensitively for the normal mass hierarchy and the uncertainties of R^ due 
to LP a , T a , and r\ a are small. Therefore, in the case of normal mass hierarchy, we can restrict #13 to a small range 
of values by measuring R\ . Similar plots are made for in Fig. 8b. Now for the normal mass hierarchy, R^ is 
independent of #13. For the inverted mass hierarchy R^ varies with 613 and is smaller than 2.4. Hence we may also 
be able to identify the mass hierarchy by However, as shown in the plot, in the case of inverted mass hierarchy 
the large uncertainties of i?^ due to T Q , and rj a make it difficult to constrain 613 by R^ . 

The Daya Bay experiment has the sensitivity limit of 0.01 for sin 2 2#i3, i.e., 3° for #13. However, if a SN explosion 
takes place during the operation of the experiment, it is possible to reach a much smaller value of 613 with the benefit 
of determining the mass hierarchy by measuring R^ for the neutrino-carbon reactions. 

Now let us consider the heavy water detector, for an example, used in the previous SNO experiment (l9j . The 
detector mass was made of 1000 tons of heavy water located in the depth of ft. = 6000 m.w.e. 4 [l~9| [5Cj . We consider only 

the neutrino-deuterium reactions. The total number of deuterium is iV^ = 6.02 x 10 31 . Scaling the experimentally 
measured energy values from the decay of the muon at rest to the e nerg y scale of SN neutrinos, one can obtain the 
effective cross sections of neutrino-deuterium reactions as follows [4!| |5l| : 

(<r(d(v e ,e-)pp)) = (3.35T 2 ; 31 - 3.70) x 10- 43 cm 2 , 

{a{d{v ei e+)nn)) = (3.05T 2 ; 08 - 7.82) x 10~ 43 cm 2 , (38) 
for the charged-current capture, and 

(a(d(v a ,v a )np)) = (1.63T 2 a 26 - 2.78) x 10 " 43 cm 2 , 

(a(d(D a , 9 a )np)) = (2.03T 2 ; 05 - 3.76) x 10 - 43 cm 2 , (39) 

for the neutral-current capture, where a = e, [A, r. 

In Figs. 9a and 9b we plot Rf and i?^ of the neutrino-deuterium reactions in analogy with Figs. 8a and 8b. Similar 
to the neutrino-carbon case given in Fig. 8a, Fig. 9a for R^ can be used to determine the neutrino mass hierarchy. 
But it is difficult to pin down the value of 6*13. Figure 9b for R± is similar to Fig. 7 and the discussion can be applied 
to the present situation, with the allowed value of R^ between 1.9 to 2.6 in the case of inverted mass hierarchy. 

In the above discussions, we have applied our method to three reactions: inverse beta decay, neutrino-carbon, and 
neutrino-deuterium. A summary is given in Table II. From this table, it can be seen that the mass hierarchy may be 
identified by measuring in the neutrino-carbon reactions and the neutrino-deuterium reactions, and by measuring 
R^ in the inverse beta decay, the neutrino-carbon reactions, and the neutrino-deuterium reactions. In addition, in 
the case of normal mass hierarchy, Rf measured in the channel of neutrino-carbon reactions can provide the most 
accurate information about 613. Similarly, the present method can also be applied to other neutrino experiments, 
such as KamLAND (52|, MinBooNE (53[, Borexino (HI, and Double-Chooz [53 in the channel of neutrino-carbon 



several authors [13 1 fill 



3 We have checked that when Rf^ and /Lj take other values, the qualitative results do not change. 

4 m.w.e. means meter-water-equivalent and 1 m of rock is about 2.7 m of water. 
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reactions. In the future, if a SN explosion takes place within the cosmic distance considered here, it is possible to 
identify the mass hierarchy and reach a much smaller value of #13 than current neutrino experiments through the 
measurement of in the channel of neutrino-carbon reactions. 



IV. CHECK NEUTRINO MASSES UNDER COSMIC SETTING 



All of our understanding of neutrino properties has been derived in the terrestrial environment, with input from 
solar neutrinos, and our study of SN neutrinos has been made under the reasonable assumption that the terrestrial 
properties can be applied unaltered to the extreme conditions of the cosmos, such as those in the SN, e.g., high 
temperature, high density, etc.. If an observation of SN neutrinos deviates from the prediction, it will be argued that 
our picture of SN requires modification. However, it would be prudent to check if some aspects of the properties of 
neutrinos may have been modified under the unusual cosmic conditions. In this section we propose a measurement, 
the observation of which, can only be attributed to some change of the neutrino property, albeit a very specific one. 

Since neutrinos, at least some of the species, have finite masses, their transmitting velocities are smaller than that 
of light. Then, neutrinos from a SN explosion, while traveling through a large distance to the Earth, can experience 
a delay in the arrival time, which depends on their energies and masses and may be observable [|[ [5(| [56| . 

Suppose a neutrino from a SN, of mass m (in eV) and energy E (in MeV), travels over a distance D (in lOkpc) to 
the Earth, it will experience an energy- dependent time delay of At (in second) relative to a massless neutrino [8[ : 

At(s) - 0.5143 [^^j D(10kpc), (40) 

where only the lowest order term in the neutrino mass expansion has been kept. In Ref. @, a method of measuring 
the neutrino mass is proposed, using the detection of the event number of the delayed SN neutrinos. In the papers, 
a simple model of the luminosities of SN neutrinos was assumed, in the absence of various effects considered in 
recent SN neutrino studies. The neutrino luminosities were assumed to have a common step-function cutoff at a 
fixed time. This form of luminosities was derived from the argument that the neutrino gravitational rcdshifts become 
severe enough only near the sharp cutoff time tsn, where tsH can be determined from the observed time profile 
of SN neutrinos. Therefore, the luminosities, which are constant before tsH and vanishes abruptly afterward, take 
the form L(t) = LP a 6(tBH — t). In the present work, we take into account the full complexity of the neutrino flavor 
transformation due to the effects of the shock, MSW, collective, and the Earth matter. Also, we consider the time- 
dependence of the SN neutrino luminosities given in Eq. (TTJ) and the uncertainties in the SN neutrino flux parameters 
given in Eqs. ([3]) and ([5]). Since the time delay At depends on the neutrino mass m, we can express the event number 
of the delayed SN neutrinos as a function of m. Therefore, information about m can be obtained by detecting the 
event number of the delayed SN neutrinos. 

From Ref. (57j there are the following neutrino mass bounds 5 

m Pe < 2eV, m Vll < 0.19MeV, m VT < 18.2MeV. (41) 

For the SN electron antineutrino the averaged energy is (En e ) ~ 16MeV, which together with D = lOfcpc, gives 
Atp e < 0.01s from Eq. (|40"|) . This time delay of the v e is much smaller than the time scale of the SN explosion of 
about 10s. And it is also too short to be be observable. Similarly, for the electron neutrino (E Ve ) ~ HMeV, we find 
At Ve < 0.02s, which is again very short and much smaller than the time scale of SN explosion. Therefore, unless 
their effective mass is somehow significantly enhanced, v e and v e can be considered as massless for the time delay 
observation. For the muon and tau neutrinos, according to the oscillation experiments their effective masses in the 
terrestrial setting will also be bound by 2eV because of the much smaller values of Am^ and Am^, 6 their time delay 
will be very small too. In the following, we will tentatively ignore these terrestrial bounds to form the basis of our 
proposal for a test of neutrino properties under the cosmic setting of extreme conditions. 



Note that these bounds are obtained from direct measurements of neutrino masses and they do not take into account neutrino oscillation 
experiments. The latter will relate neutrino effective masses of different flavors in a very tight range. 

The effective masses of the flavor neutrinos, m a , a = e, fi, and r, to those of the regular neutrino masses, mj, j = f , 2, and 3, by 

m a = \U a j\ 2 m,j, 
3=1,2,3 

where U a j is the neutrino mixing matrix. 
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After the SN explosion the emission of different flavors of neutrinos lasts for a period of about 10s. The arrival 
to the Earth of those flavors of neutrinos and antineutrinos whose effective masses are not affected will last about 
the same period of time. For other flavors which are affected there may be a delay in the arrival time, depending 
on the values of their effective masses. We look for measurements which are sensitive to both charge and neutral 
current events, such as the detector types of the Daya Bay [H, LVD 58 1, and SNO [59[ experiments, which involve the 



neutrino-carbon and neutrino-deuterium reactions. The experimental information that these detectors can potentially 
provide include the event numbers of the electron antineutrino from charged-current reactions and the total event 
number of all the three flavors of neutrinos from the neutral-current reactions. We shall assume that the electron 
neutrino mass is not affected so as to avoid any possibility of influencing the existing result of nucleosynthesis. But 
the effective masses of the muon and tau neutrinos might be affected. Then, there might be a time delay At x between 
the charged-current events and some of the neutral-current events if the muon and tau flavors of SN neutrinos appear 
to have sufficiently sizable effective masses. Denoting the effective mass of v x (i/ M , 9^, v T , or D T ) by m„ x , we have 

At x = (t) NC -(t) cc ^Q.5U3(^J D, (42) 

where we have dropped the units for At x , m x , E, and D as indicated in Eq. (j40|) . < t >nc is the average time of SN 
neutrinos of the neutral-current reactions to arrive at the detector and < t >cc that of the charged-current reactions. 
A nonvanishing At x signalizes a significant larger effective mass for muon or tau neutrino, independent of parameters 
of the SN and neutrinos. 

To gain a sense of how At x may look like and how sensitive it is to the effective mass of the muon and tau neutrinos, 
we have to do some model calculation for quantities which can be measured in a SN neutrino experiment. We found 
that such a quantity is available as described below. We use the charged-current events which are from the electron 
antineutrino to monitor the events that are not delayed, and the neutral-current events which involve all flavors to 
measure the possible time delay. We divide the total neutral-current events into two groups: N^PJf and N^Jf\ 
where N^Jf^ is the total neutral-current events in accompanying with the charged-current events and are therefore 

not delayed, and N^fjf^ are all the total neutral-current events including both delayed and not delayed. We define 
the event number of the delayed v x plus v x from the neutral-current events as 

The quantity to be measured is the new ratio, 

N™ (delay) Agg) _ 

2 N NC(0) M NC(0) ■ 

ls total ls total 

This ratio is very insensitive to the values of Am\ x and Ato 2 ^. 7 It should also be remarked that because of the ratio, 
the effect of the various uncertainties involved, such as the neutrino and SN parameters, will be reduced in R2. Using 
the charged current events to monitor the undclaycd neutral current events will bypass the question of the uncertainty 
in the total emission time of SN neutrinos. 

Since we can not tag the flavor of neutral-current events the information on the masses of the muon and tau 
neutrinos are correlated. For the purpose of illustration we assume and v T to have the same effective mass. We 
keep all other features of neutrinos and the SN unchanged. For the neutrino-carbon reactions at Daya Bay, we plot 
in Fig. 10 two sets of curves which describe R2 as a function of m Ux , each for a fixed value of #13. In Fig. 10(a) 
for = 0, it can be seen that the curves are insensitive to L° jlP v and the mass hierarchy. For m„ x < 200eU, 
i?2 is sensitive to m^. A measurement of R2 can restrict m x in a distinctive range of values determined by the 
uncertainties in the SN neutrino temperature and chemical potential. In Fig. 10(b) for 8 13 = 3°, R2 is again shown 
to be sensitive to m x for m Vsc < 200eV. The two mass hierarchies overlap for small m x but fall into well separated 
regions for m. x > 200^. 

In Figs. 11 we plot R2 for the neutrino-deuterium reactions in analogy to Figs. 11 of the neutrino-carbon reactions. 
Figure 11(a) is very similar to Fig. 10(a). Figure 11(b) is also similar to Fig. 10(b) for m Ux smaller than 150eV, but 
there are some differences at larger m x . In this latter case, unlike Fig. 10(b), Fig. 11(b) shows that the two mass 
hierarchies can no long be separated even for m„ x greater than 200eU. 



7 Our numerical simulation shows that R2 is insensitive to variations of Arrijj and Arrtgj in a range of 10 and 8 orders of magnitude 
respectively. 
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The method can also be applied to KamLAMD, MinBooNE, Borcxino, and Doublc-Chooz in the channel of neutrino- 
carbon reactions, and the results are similar to those given for the Daya Bay experiment. 



V. DISCUSSION AND SUMMARY 



In this paper, we considered the SN shock effects, the MSW effects, the collective effects, and the Earth matter 
effects in the detection of type II SN neutrinos. Also, we considered the uncertainties in the neutrino luminosities L° a , 
their temperatures T Q , and their pinching parameters ry Q in the calculation of different flavors of SN neutrinos. We 
found that quantities suitably defined in terms of event numbers of different flavors of SN neutrinos are sensitive to 
the neutrino mass hierarchy, the mixing angle #13, and neutrino masses. Therefore, it is possible to identify the mass 
hierarchy, acquire information about $13, and look into certain neutrino mass patterns of SN neutrinos. 

Firstly, we defined the ratio R^~^ as the event number of v e (y e ) for the charged-current reaction detected after 

one second over that detected in the first one second, and discussed the dependence of RX on 613 for different 
mass hierarchies. With this, we may identify the mass hierarchy and obtain information about #13. This method can 
be applied to several detector types relying on specific physical processes, including the inverse beta decay (Fig. 7), 
neutrino-carbon reactions (Figs. 8(a) and 8(b)), and neutrino-deuterium reactions (Figs. 9(a) and 9(b)). From our 
numerical calculations, the mass hierarchy may be identified by measuring R^ in the neutrino-carbon reactions and the 
neutrino-deuterium reactions, in the inverse beta decay, the neutrino-carbon reaction, and the neutrino-deuterium 
reactions. In addition, in the case of normal mass hierarchy, R^ measured in the neutrino-carbon reactions provides 
the most restrictive information about #13. Therefore, we can both identify the mass hierarchy and measure #13 from 
R^ in the channel of neutrino-carbon reactions at several experiments such as Daya Bay, KamLAND, MinBooNE, 
Borcxino, and Double-Chooz. 

Secondly, we defined another ratio, R2, for the neutral-current reactions, as the event number of the delayed neu- 
trinos over the total event number, with the charge current events as a monitor of undclaycd events. The significance 
of i?2 is that it offers a possibility to detect possible difference between the electron neutrino and its other flavors 
counterparts, in particular, any significant difference in their effective masses under the cosmic setting of extreme 
conditions. With the plot of Ri vs m„ x , we can study m Vx vs to„ c . This method is applied to both neutrino-carbon 
and neutrino-deuterium reactions and the result is given respectively in Figs. 10 and 11. What we have presented for 
i?2 is a case study with rather restricted input. In case there is an indication of non-vanishing Ri a detailed refined 
calculation should be made. 

In conclusion, if a SN explosion takes place within the cosmic distance considered here, it is possible to identify 
the mass hierarchy and obtain information on #13. With suitable set up of the detector, a study of some unusual 
properties of neutrinos under the cosmic setting is possible. A summary of the detector types and their capabilities 
in association with the studies presented in this article is given in Tabic II. 
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Figure Captions 

Fig. 1 The density of the SN as a function of the radius r at different times. 

Fig. 2 The crossing probability at the high resonance region Pg as a function of the neutrino energy E for two very 
small neutrino mixing angles #13 at t = 4s. 

Fig. 3 The crossing probability at the high resonance region Ph as a function of time t for three different neutrino 
energies at sin 2 #13 = 10~ 4 . 

Fig. 4 The crossing probability at the high resonance region Pjj as a function of the neutrino mixing angle 6* 13 : (a) for 
three different times at E = HMeV. The solid, dashed, and dotted curves are for t = 2s, 4s, 6s, respectively; (b) for 
three different neutrino energies at t = 4s. The solid, dashed, and dotted curves are for E = HMeV, 16MeV, 25MeV, 
respectively. 

Fig. 5 Illustration of the path of the SN neutrino reaching the detector in the Earth. D is the location of the detector, 
6 is the incident angle of the neutrino, O is the center of the Earth, L is the distance the neutrino travels through the 
Earth matter, and x is the instantaneous distance of the neutrino to the center of the Earth. 

Fig. 6 The luminosity decay time t as a function of the mixing angle #13 in the channel of neutrino-carbon reactions 
at Daya Bay when = 2.4 and Li J Li = 1. The "max" ("min") corresponds to the maximum (minimum) values 
of T a and r/ a . 

Fig. 7 The ratio of the event number of D e for the charged-current reaction detected after Is over that detected 
before Is, as a function of the mixing angle #13 in the channel of inverse beta decay at Daya Bay. The solid curves 
correspond to L® /L® =1/2, the dashed curves Li / Li =1, and the dot-dashed curves Li /Li = 2. The "max" 
("min") corresponds to the maximum (minimum) values of T a and r\ a . 

Fig. 8 (a) Similar to Fig. 7 but the electron neutrinos v e are observed in the channel of neutrino-carbon reactions 
at Daya Bay; (b) similar to Fig. 7 but the anti-electron neutrinos D e are observed in the channel of neutrino-carbon 
reactions at Daya Bay. 

Fig. 9 (a) Similar to Fig. 7 but the electron neutrinos v e are observed in the channel of neutrino-deuterium reactions 
at SNO; (b) similar to Fig. 7 but the anti-electron neutrinos v e are observed in the channel of neutrino-deuterium 
reactions at SNO. 

Fig. 10 For the neutral-current reactions R2 in the channel of neutrino-carbon reactions at Daya Bay, the ratio of 
the event number of delayed neutrinos over the total event number, as a function of m Vi at two values of #13, each 
has the two mass hierarchies with Li /L® , T Q , and r) a taking their limiting values. The solid curves correspond 
to L°jL° x = 1/2, the dashed curves i^/i^ = 1, and the dot-dashed curves L a v JL a Vx = 2. The "max" ("min") 
corresponds to the maximum (minimum) values of T a and rj a . Figure (a) shows that R2 is nearly independent of 
Li I Li and the mass hierarchy. 

Fig. 11 Similar to Fig. 10 but in the channel of neutrino-deuterium reactions at SNO. 
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TABLE I: Summary of neutrino flavor conversions due to the neutrino shock effects and the MSW effects in various density 
regions. 



t < Is 


t > Is 




Pb > p- 


p 


- > Pb 




resonance in 
region 


flavor conversion 
involved 


Resonance in 
region 


flavor conversion 
involved 


Ph3 


Prea > P+ 


Phi, Pa 


Prea > P+ 


Phi, Pa 




P+ > Prea > Pb 


Phi, Ph2, Pa 


P+ > Prea > P- 


Phi, Ph2, P s 




prea = Pb 


Phi, Pa 


P- > Prea > Pb 


Phi, Phi, Pa, Ph 3 




Pb > Prea > p- 


Pa 


prea = Pb 


Phi, Pa Phs 




P- > Prea 


Pa, P H3 


Pb > prea 


Pa, PH3 



TABLE II: Summary of measurement of the neutrino parameters by SN neutrinos in three reactions: the inverse beta decay, 
neutrino-carbon reactions, neutrino-deuterium reactions, "y/" and "x" represents respectively the information can and can 
not be obtained. 
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Fig. 7 
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